Effect of dimensionality on the charge-density-wave in few-layers 2H-NbSe2 
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We investigate the charge density wave (CDW) instability in single and double layers, as well as 
in the bulk 2H-NbSe2. We demonstrate that the density functional theory correctly describes the 
metallic CDW state in the bulk 2H-NbSe2- We predict that both mono- and bilayer NbSe2 undergo a 
CDW instability. However, while in the bulk the instability occurs at a momentum qcow ~ §rM, 
in free-standing layers it occurs at <\cdw ~ |rM. Furthermore, while in the bulk the CDW leads 
to a metallic state, in a monolayer the ground state becomes semimetallic, in agreement with recent 
experimental data. We elucidate the key role that an enhancement of the electron-phonon matrix 
element at q w (\cdw plays in forming the CDW ground state. 
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Charge density wave (CDW) is one of the most 
common and most intriguing phenomena in solid state 
physics The concept originated in the seminal paper 
by PeierlsQ pointing out a divergence in the ID response 
functions at a wave vector equal to twice the Fermi vec- 
tor. Later this concept was generalized onto the 2D and 
3D system with "nesting", that is, quasi- ID portions of 
the Fermi surface. On the other hand, the fact that the 
Peierls instability is logarithmic, and therefore very frag- 
ile, has led to the question of whether the actual CDW 
observed in quasi-2D materials are indeed manifestations 
of the Peierls instability. Q 

In the last decades the role of nesting in various 
structural and magnetic instabilities has been wildly de- 
bated, particularly with respect to the most venerable 
CDW materials, NbSe2 and related dichalcogenides. It 
has been pointed out that the bare susceptibility does 
not have any sharp peak at the CDW wave vector, 
<Icdw_ = (2tt/3, 2tt/3), but at best a broad and shallow 
peakjj, 0, H[ in the real part, while the imaginary part 
(which directly reflects the Fermi surface nesting) does 
not peak at qczjw at all[l, 

This suggests that momentum dependence of the 
electron-phonon interaction plays a crucial role in driv- 
ing CDW instabilities]!, @] . An indirect confirmation 
comes from neutron (8j| and X-ray scattering 0] phonon- 
dispersion data where a marked softening was detected 
in only one of the low energy modes close to c\cdw- 

A very clean test for the described hypothesis are 
purely 2D systems — mono- and bilayers of the same 
compounds. Lacking k z dispersion, they should be much 
more sensitive to the Fermi surface nesting than their 3D 
counterparts. 2H-NbSe2 is particularly promising as its 
2D structures have been synthesized, and found to show 
a curious semimetallic behavior [Io| . Interestingly the 
conductivity of a NbSe2 layer in a field-effect transistor 
as a function of the gate voltage (V g ), which controls the 
doping, is very different from that of good metallic mono- 



layers. While typical metals have a conductivity rather 
weakly dependent on doping, NbSe2 monolayers show a 
strong monotonic dependence on V g . Such a large varia- 
tion (more then a factor of two in the -70 to +70 V V g - 
range) is characteristic of semiconductors, semimetals as 
graphene, or metals with a pseudogap at the Fermi level. 
The first principles calculations predict a good metal for 
a monolayer in the undistorted bulk structure [ll[ , in dis- 
agreement with experiments, which suggest a possibility 
of a semimetallic band structure. 

In this Letter we study the CDW instability in the bulk 
2H-NbSe2, in the bilayer and in the monolayer NbSe2, 
both by calculating the phonon spectra in the high- 
symmetry phase and by full structural optimization[l2|. 
We find that all three studied structures are unstable 
against a CDW formation with, however, different pat- 
terns, despite the extreme similarity of the Fermi sur- 
faces. We then demonstrate that the crucial difference 
comes from the q-dependent electron-phonon matrix ele- 
ments, rather than directly from the Fermi surface (FS). 
Finally we show that the CDW distorted monolayer is 
semimetallic and the large variation of its conductivity 
as a function of V g detected in experiments is a manifes- 
tation of the CDW. 

We simulate a single layer or a bilayer, as usually, by 
introducing a thick vacuum layer (of the order of 10 A). 
We then performed the geometrical optimization in one 
unit cell, without a CDW. We found that the distance 
between the Nb and Se planes barely changes: from 1.67 
A in the bulk to 1.69 A in the bilayer (same in the 
monolayer) . On the contrary, the Nb-Nb interplanar dis- 
tance in the bilayer is significantly enlarged compared to 
the bulk (6.27 A vs. 6.78 A). 

The pseudopotential electronic structure of the NbSe2 
bilayers and monolayers is compared with that of the 
bulk NbSe2 in Fig. [T] (LAPW calculations yield very 
similar results). For the bulk CclSC, ctS known from previ- 
ous calculations [H, HH , three bands cross the Fermi level 
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FIG. 1: Electronic structure of the undistorted bulk, bilayer 
and monolayer 2H-NbSe2. 




FIG. 2: (color online) Left: Real part of the bare static elec- 
tronic susceptibility (in arbitrary units) with constant ma- 
trix elements, same as the real part of the phonon self-energy 
II . iAt (q,w = 0), Eq. (fl), calculated using the LAPW band 
structure of a single layer NbSe2- Center: The corresponding 
FS. Four Brillouin zones are shown, V points are at the grid 
nodes. Right: The same, for the resulting 4x1 CDW. 



(ef). One is the antibonding Se p z bands, which is highly 
dispersive along the z-axis and does not cross et in the 
A-L-H plane. The other two are formed by Nb d-states. 
These two bands are degenerate in the A-L-H plane in 
the scalar relativistic calculations (spin-orbit interaction 
removes the degeneracy) while the degeneracy is lifted for 
other values of k z [B[ by the interlayer interaction. As a 
result, while one of the FSs originated from Nb d-bands 
is two dimensional, the other acquires considerable k z 
dispersion, so that a model neglecting hopping between 
the layers does not apply [l6j| . 

In a bilayer, the Se-Se interplanar interaction is re- 
duced, while in the monolayer it is absent. As a result, 
the Se band sinks below the Fermi level, and only Nb the 
derived Fermi surfaces are present, and they are strictly 
2D, that is, seemingly more prone to the nesting effects. 
The calculated FS for a single layer (Fig. [2]) consists of a 
rounded hexagon at T and two rounded triangles at K, 
just as the Nb FS pockets in the bulk 2H-NbSe 2 d . 

Next, we have computed the phonon dispersions for 



the bulk and single-layer NbSe2 . In Fig. [3] we plot the 
phonon dispersions in the monolayer (bulk) calculated 
using a 40 x 40 (20 x 20 x 6 ) k-points grid and Fermi 
temperature r = 68 meV (t = 270 meV) with Hermite- 
Gaussian distributions (See Ref. [l7[ for more details.). 
In the bulk, the highest energy acoustic mode is unstable 
at approximately qcDW — 2/3rM, in agreement with 
the experiment. These calculations implicitly include the 
renormalization of the bare phonon frequencies due to the 
electron-phonon interaction, as described by the real part 
of the phonon self-energy of a phonon mode /x, namely 
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where £ki is the one-electron energy, k+q , is the 
electron-phonon matrix element and is the Fermi 
function. If one neglects the k dependence of the matrix 
elements in Eq. [TJ this expression becomes proportional 
to the real part of the unrenormalized static electronic 
susceptibility, x'oi'i)'- as shown in Ref. [H, in NbSe2 Xoi'i) 
has a broad maximum at (\cdw = 2/3rM, even though 
the imaginary part (which directly reflects the FS nest- 
ing) does not. 

Thus, DFT correctly describes the CDW instability 
in the bulk NbSe2. With this in mind, we move to a 
NbSe2 monolayer and observe that (a) now the most 
unstable phonon mode along TM appears near (l/2)rM 
= (0.57r/a,0) and (b) the instability expands over a sub- 
stantially larger k-point region, as phonon are unstable 
also along MK. We also looked at the bare susceptibility 
Xo and found a strip- like region of enhanced values, ex- 
tending from 0.4rM to 0.75FM, with a width of approxi- 
mately 0.15 of the TM vector (0.157r/a). This is in agree- 
ment with the instabilities found in the linear-response 
phonon calculations along TM and MK. 

To understand which electronic states are involved 
in the CDW formation we have repeated the calcula- 
tions with an increased electronic temperature r. We 
found that the phonon frequency of the highest acoustic 
mode strongly depends on r for r > 0.3eV. The insta- 
bility occurs for r rj 0.3 eV, and for smaller tempera- 
tures the (imaginary) phonon frequency of the unstable 
mode is essentially constant. This means that electronic 
states within a 0.3 eV window around e/ are responsi- 
ble for the CDW instability, which is not a Fermi sur- 
face effect [1, HH, EH, EEl, but involves states in the full 
Nb-derived d-band (which is only 1 eV wide, see Fig. 
|4}. As discussed above, x'o does have a maximum at 
QCDWi but this maximum is relatively weak and broad, 
and, most importantly, weakly dependent on r. Thus, 
the CDW is mainly the result of an enhancement of the 
electron-phonon matrix element close to <\cdw ■ The self- 
consistent screening of the electron-phonon matrix ele- 
ment as a function of r is crucial to describe the CDW. 

The next test is to take a superstructure suggested by 
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FIG. 3: (color online) Phonon dispersion of bulk (top) and 
monolayer (bottom) 2H-NbSe2. Each point corresponds to a 
linear response calculations. The spline connecting the points 
is guide to the eye. 



the calculated CDW vector [(tt/3,0) or (tt/4,0)] and to 
optimize the cell dimensions and the atomic positions. 
Each time we started from slightly randomized atomic 
positions and minimized the total energy. We found that 
in the bulk the 4x1 and 3x1 supercells both converge to 
a lower energy than the undistorted cell, consistent with 
the calculated phonon spectra in Fig. 02 the 3x1 being 
lower in energy (with the energy gain of 0.1 mRyd/Nb 
with respect to the undistorted structure). The other 
supercells converged to the undistorted structure. The 
distortion of the 3x1 supercell is consistent with the £i 
honon pattern detected in inelastic neutron scattering 
, involving not only an in-plane deformation of the lay- 
ered structure, but also out-of-plane displacements of the 
Se atoms and a charge transfer between the two layers. 

For a monolayer, the situation is different. The lowest 
energy supercell is now 4x1 with a 1.19 mRyd/Nb energy 
gain with respect to the undistorted cell, followed by 3x 1. 
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FIG. 4: (color online): Pseudopotentials electronic density 
of states for undistorted (black-dashed) and distorted (red) 
bulk, bilayer and monolayer 2H-NbSe2. 



We also checked the 2x1 and 2x2 supercells, which 
converged to the undistorted structure, and the 3x3 one, 
which converged to the same structure as 3 x 1. This is 
again consistent with the calculated phonon spectra. In 
the ground state Nb atoms are trimerized, qualitatively 
different from the bulk CDW. Finally, the bilayer behaves 
just like a monolayer, with a 4 x 1 ground state similar 
to that found for a monolayer. 

Having obtained the distorted structure both for the 
bulk and for a monolayer, we investigated their electronic 
properties. We found (Fig. [4]) that in the bulk case the 
density of states (DOS) at the Fermi level in the dis- 
torted 3x1 supercell is not much reduced in the CDW, 
contrary to what one would expect in the Peierls model, 
and reconfirming that the energy gain in this case is not 
coming from the states in the immediate vicinity of the 
Fermi levelQ but from the occupied Nb bands. In the 
monolayer 4x1 supercell the Peierls model is again vi- 
olated due to the important role of the electron-phonon 
matrix element, however the FS is strongly gapped (see 
Fig. [5] right) and displays a clear semimetallic behavior 
in qualitative agreement with experimental data (loj . 

More insight on the conducting nature of the mono- 
layer is obtained from the conductivity per layer, namely 



2D 



(E) 



3n 2D N, 



ki 



vl x T(k)S(e ki - E). (2) 



where Q2D is the area of the two dimensional unit cell of 
a NbSe2 layer, e is the electron charge, Nk the number of 
k-points used in the calculation, Vki the electron velocity, 
E the energy. In general a^ c \E) — a^P{E) x L, where 
L = 6.275A is the layer thickness. Assuming, r(k) w 
r( e/ ), we have that a™ (E) = T{e } )u 2 p ab i 'Air . 
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FIG. 5: (color online): LAPW plasma frequency and con- 
ductivity as a function of the gate- induced charge for a single- 
layer NbSe2- The constant relaxation time approximation is 
adopted. Experimental data [10] are plotted assuming the 
charge neutrality at ~ 30 x 10 12 cm -2 . The gate voltage 
for the standard configuration of a 300 nm SiC>2 field-effect- 
transistor can be obtained from the x axis dividing 7.2 x 10 -2 
cm -2 V -1 . The charge is obtained from the integrated DOS 



The results of the ui^ ab calculation using the optics 
code of the WIEN2k package are plotted in Fig. [5] for 
the monolayer NbSe2 as a function of the gate-induced 
charge as obtained from the calculated DOS (Fig. Q|. 
For the undistorted structures, uj p ,ab is essentially con- 
stant, in disagreement with the experiment. On the con- 
trary, in the CDW state the monolayer plasma frequency 
shows substantial variation near ep, due of its semimetal- 
lic structure. A pseudogap occurs just above e/, so that 
the conductivity grows with the gate voltage. In the 
constant relaxation-time approximation, and assuming 
an electron self-doping of w 30 x 10 12 cm~ 2 we obtain 
reasonable agreement with experimental conductivity as 
function of the gate volt age [10] (see Fig. O. For a quanti- 
tative comparison with experiments, however, one would 
need to know the energy and momentum dependence of 
the relaxation time. 

To summarize, we have shown that the density func- 
tional theory accurately describes the CDW instability 
in the bulk 2H-NbSe2, and predicts a similar instabil- 
ity to occur in mono and bilayer. However, while in the 
bulk the CDW occurs at qcDW ~ (§,0,0)^, in the 
monolayer and bilayer the CDW vector is w (^,0) — ■ 
In all these systems the CDW is driven by an enhance- 
ment of the electron-phonon coupling at q ~ qc dw ■ Our 
work solves the long standing controversy on the origin 
of CDW in 2H-NbSe2 and in transition metal dichalco- 
genides. 

Unlike the bulk case, when the system remains a good 



metal in the CDW state, in a monolayer CDW produces 
a semimetallic state. This leads to a large variations of 
the conductivity as a function of the gate voltage, con- 
trary to the undistorted structure and in agreement with 
the experiment data(Io|. Our calculation shows that a 
proper description of the CDW instability in reduced 
dimension is mandatory to interpret conduction data 
for low-dimensional transition-metal-dichalcogenides and 
opens the way to theoretical understanding of CDW- 
based field-cffcct-dcvices. 



Calculations were performed at the IDRIS supercom- 
puting center (project 081202). 
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